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The recent ALEPH measurements of the inclusive Cabibbo-suppressed decay width of the r and several moments 
of its invariant mass distribution are used to determine the value of the strange quark mass. We obtain, in the 
MS scheme, m 3 (M 2 ) = (119 ± 24) MeV, which corresponds to ro s (l GeV 2 ) = (164 ± 33) MeV, m s (4GeV 2 ) = 
(114 ±23) MeV. 



1. INTRODUCTION 

The precise numerical value of the strange 
quark mass is a controversial issue, with impor- 
tant implications for low-energy phenomenology. 
The Particle Data Group M quotes a rather wide 
range of m s values, reflecting the large uncertain- 
ties in the present determinations of this param- 
eter from QCD Sum Rules and Lattice calcula- 
tions. 

The high precision data on tau decays col- 
lected at LEP and CESR provide a very powerful 
tool to analyse strange quark mass effects in a 
cleaner environment. The QCD analysis of the 
inclusive tau decay width, 

T [t~ — > v T + hadrons (7)] 
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has already made possible |3[| an accurate mea- 
surement of the strong coupling constant at the 
r mass scale, a s (M%), which complements and 
competes in accuracy with the high precision 
measurements of a s (M§) performed at LEP. 
More recently, detailed experimental studies of 
the Cabibbo-suppressed width of the r have 
started to become available |i|,|5| , allowing to per- 
form a systematic investigation of the corrections 
induced by the strange quark mass in the r decay 
width §-§. 

What makes a m s determination from r data 
very interesting is that the hadronic input does 
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not depend on any extra hypothesis; it is a purely 
experimental issue, which accuracy can be sys- 
tematically improved. The major part of the un- 
certainty will eventually come from the theoret- 
ical side. However, owing to its inclusive char- 
acter, the total Cabibbo-suppressed tau decay 
width can be rigorously analyzed within QCD, 
using the Operator Product Expansion (OPE). 
Therefore, the theoretical input is in principle un- 
der control and the associated uncertainties can 
be quantified. 

2. THEORETICAL FRAMEWORK 

The theoretical analysis of the inclusive 
hadronic tau decay width [p[-|T2]| involves the two- 
point correlation functions 



i d'xe^ (0\T(j$(x)J$(0)i) |0) 



for the vector, J- 1 - = V^(x) = qfl^qi, and axial- 
vector, j£ = A(j(x) = ^7^75^, colour-singlet 
quark currents — u,d,s). These correlators 
have the Lorentz decompositions 



ij,V/A 
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where the superscript in the transverse and lon- 
gitudinal components denotes the corresponding 
angular momentum J = 1 (T) and J = (L) in 
the hadronic rest frame. 
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Table 1 

Explicit values of the relevant kinematical kernels. 



The semi-hadronic decay rate of the r lepton, 
can be expressed as an integral of the spectral 
functions Imil T (s) and lmIl L (s) over the invari- 
ant mass s of the final-state hadrons as follows: 



Mi 



R t = 12tt — 1 
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MS 
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Moreover, according to the quantum numbers 
content of the two-point function correlators 



U J (s) 



\Vu d \ 2 [ni ud (s) + KuM 



IK S | 2 [n J v, us (s) 



we can decompose R T into 

R T = R T y + R T ^A + Rt,S , 



(4) 



(5) 



where R T ,v and R t ,a correspond to the first two 
terms in Eq. (^J), while R T ,s contains the remain- 
ing Cabibbo-suppressed contributions. 

The measurement of the invariant mass distri- 
bution of the final hadrons provides additional 
information on the QCD dynamics, through the 
moments [E2| 



which include R T = R®° as a particular case. 

Exploiting the analytic properties of n J (s), we 
can express these moments as contour integrals 
in the complex s-plane running counter-clockwise 
around the circle \s\ — M 2 : 



R k } = -m 



^{?,Ti l +T {x)D L + T {M*x) 



We have used integration by parts to rewrite Rf 
in terms of the logarithmic derivatives 
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d l ( s ) = *—\ 8 n L (s)] , 

Ml ds L 
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which satisfy homogeneous renormalization group 
equations. All kinematical factors have been ab- 
sorbed into the kernels T^ +T (x) and J 7 ^ l (x). Ta- 
ble [l] shows the explicit form of these kernels for 
the moments which we are going to analyze in the 
following sections. 

For large enough — s, the contributions to 
D J (s) can be organized with the OPE in a se- 
ries of local gauge-invariant scalar operators of 
increasing dimension D — 2n, times the appro- 
priate inverse powers of — s. This expansion is 
expected to be well behaved along the complex 
contour \s\ = M 2 , except in the crossing point 
with the positive real axis |13|. As shown in Ta- 
ble |, the region near the physical cut is strongly 
suppressed by a zero of order 3 + k at s = M 2 . 
Therefore, the uncertainties associated with the 
use of the OPE near the time-like axis are very 
small. Inserting this series in (Q) and evaluating 
the contour integral, one can rewrite Rlf as an 
expansion in inverse powers of M 2 Q , 

Rf = 3 [\V ud \ 2 + \V US \ 2 ] S E w (l + <W + *«(°> 
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where sin 2 6n = 



\V US \ 2 l\\V ud \ 2 



we have pulled out the electroweak corrections 



• (7) Sew = 1-0194 (lj and S' EW ~ 0.0010 | 



The dimension-zero contribution S kl ^ is the 
purely perturbative correction, neglecting quark 
masses, which, owing to chiral symmetry, is iden- 
tical for the vector and axial-vector parts. The 
symbols S^ D) = [5%™ + 6%$ ] ]/2 stand for 
the average of the vector and axial-vector con- 
tributions from dimension D > 2 operators; they 
contain an implicit suppression factor 1/M^jP. 

3. SU(3) BREAKING 

The separate measurement of the Cabibbo- 
allowed and Cabibbo-suppressed decay widths of 
the r [fi} allows one to pin down the SU(3) break- 
ing effect induced by the strange quark mass, 
through the differences 
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The leading contributions to SR'f are quark- 
mass corrections of dimension two HQ ; they are 
the dominant SU(3) breaking effect, generating 
the wanted sensitivity to the strange quark mass. 
The corrections of 0(m 4 ) are very tiny j7j. The 
main D = 4 contribution comes from the SU(3)- 
breaking quark condensate 

S0 4 = (0\m s ss-m d dd\Q). (11) 

Neglecting the small (9(m 4 ) terms and D > 6 
contributions, SRf can be written as H: 



SR. 



24 S 



EW 



Ml 



(1-4) Ag)(a T ) 



2 „ i s, 

^ T ' ) 



(12) 

= 0.053 ±0.002 [|16[ and a T = 

(2) 



where = wid/m 
a s (M T 2 )/7r. 

The perturbative QCD expansions A^(a 
and Qki{a T ) are known to 0{a 2 .). Moreover, the 
O(a^) contributions to A^(o T ) coming from the 
longitudinal correlator D L (s) have been also com- 
puted. Using the value of the (MS) strong cou- 
pling determined by the total hadronic r decay 
width §, a s {M?) = 0.35 ± 0.02, one gets the 
numerical results shown in Table || jl]] . 
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Table 2 

Numerical values of the relevant perturbative 
expansions for a s (M?) = 0.35 ± 0.02. 

The rather large theoretical uncertainties of 



AgV) = - A {3A^+ T K) + A^K)} , (13) 

have their origin in the bad perturbative be- 
haviour of the longitudinal contribution. The 
most important higher-order corrections can be 
resummed 0, using the renormalization group, 
but the resulting "improved" series is still rather 
badly behaved. For instance, 

A^j(O.l) = 1. 5891 + 1. 1733 + 1. 1214 + 1.2489+ ■• • 

which has 0(a 2 ) and 0(a 3 ) contributions of the 
same size. On the contrary, the J = L + T series 
converges very well: 

A£+ T (0.1) = 0.7824 + 0.2239 + 0.0831 + • • • 

Fortunately, the longitudinal contribution to 

(2) 

A kl (a T ) is parametrically suppressed by a factor 
1/3. Thus, the combined final expansion looks 
still acceptable for the first few terms: 



1 



A$ (0.1) = 0.9840 + 0.4613 + 0.3427 
+ (0.3122 - 0.000045 c^ +T ) + • • 



(14) 



Nevertheless, after the third term the series ap- 
pears to be dominated by the longitudinal contri- 
bution, and the bad perturbative behaviour be- 
comes again manifest. Taking the unknown 0(a 3 ) 
coefficient of the D L+T (s) perturbative series as 

C 3 +T ~ C 2 +T ( C 2 +T / C i +T ) ~ 323, the fourth 
term becomes 0.298; i.e. a 5% reduction only. 

Since the longitudinal series seems to reach an 
asymptotic behaviour at 0(a 3 ), the central val- 

(21 

ues of AJy (a T ) have been evaluated adding to 



the fully known 0(a 2 ) result one half of the lon- 
gitudinal 0(a 3 ) contribution. To estimate the as- 
sociated theoretical uncertainties, we have taken 
one half of the size of the last known perturba- 
tive contribution plus the variation induced by a 
change of the renormalization scale in the range 
£ G [0.75, 2] (added in quadrature). 

The SU(3)-breaking condensate SO4 could be 
extracted from the r decay data, together with 
m s , through a combined fit of different 5R* mo- 
ments. However, this is not possible with the 
actual experimental accuracy. We can estimate 
the value of SO4 using the constraints provided 
by chiral symmetry. To lowest order in Chiral 
Perturbation Theory, 5O4 is fully predicted in 
terms of the pion decay constant and the pion 
and kaon masses: SO4 ~ — / 2 (mj { — m 2 ) ~ 
-1.9 x 10~ 3 GcV 4 . Taking into account the lead- 
ing 0(p 4 ) corrections through the ratio of quark 
vacuum condensates Jl7| , ^8| 

v. S $m =0.8 ±0.2, (15) 
(0\dd\0) V ' 

one gets the improved estimate, 
S0 4 ^ -rf {v s -td) flml 

~ -(1.5 ±0.4) x lO^GeV 4 , (16) 

where we have used the known quark mass ratio 
@ m s /m = 24.4 ± 1.5. 

Strictly speaking, SO4 and v s are scale depen- 
dent. This dependence cancels with the 0(m 4 ) 
contributions and is then of 0(p 8 ) in the chiral 
expansion. The numerical effect is smaller than 
the accuracy of ( |l6| ) and has been neglected to- 
gether with the tiny 0(m 4 ) corrections. 

4. NUMERICAL ANALYSIS 

The ALEPH collaboration has measured Q the 
weighted differences SRf for five different values 
of (k,l). The experimental results are shown in 
Table g, together with the corresponding m s (M 2 ) 
values. Since the QCD counterparts to the mo- 
ments (fc, I) — (1,1) and (1,2) have theoretical un- 
certainties larger than 100%, we only use the mo- 
ments (k,l) = (0,0), (1,0), and (2,0). 

The experimental errors quoted in Table || do 
not include the present uncertainty in \V US \. To 
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Table 3 

Measured Q moments 5R^ 1 and corresponding 
to s (M 2 ) values Jt|. The first error is experimental 
and the second theoretical. 



estimate the corresponding error in m s , we take 
the following numbers from ALEPH: R^°y +A = 
3.486 ±0.015, i?°° s = 0.1610 ± 0.0066, \V ud \ = 
0.9751 ±0.0004 and \V US \ = 0.2218 ± 0.0016. This 
gives SRf = 0.394 ± 0.135 ± 0.047, where the 
second error comes from the uncertainty in \V US \ 
and translates into an additional uncertainty of 
10 MeV in the strange quark mass. We will put 
the same \V US \ uncertainty to the other two mo- 
ments, for which the ALEPH collaboration does 
not quote the separate values oiRf v+A and Rf s - 
Taking the information from the three mo- 
ments into account, we get our final result |Q: 

m s (Ml) = (119 ±12 ±18 ±10) MeV 

= (119 ±24) MeV. (17) 

The first error is experimental, the second reflects 
the QCD uncertainty and the third one is from 
the present uncertainty in | V^ s | . Since the three 
moments are highly correlated, we have taken the 
smaller individual errors as errors of the final av- 
erage. Our determination ( |l7j ) corresponds to 

m s (lGeV 2 ) = (164 ±33) MeV (18) 

and 



m s (4GeV 2 ) = (114 ±23) MeV. 



5. COMPARISON WITH ALEPH 



(19) 



The ALEPH collaboration has performed a 
phenomenological analysis of the SR 1 ^ 1 moments 
in Table |^, which results in larger m s values Q : 

r.ris f 149 + 3o cxp+ 2^ ± 6fit MeV, 

T> \ 176+^-l 2 2 4 -±8 fit ±ll J=0 MeV. 



To derive these numbers, ALEPH has used our 
published results in refs. ||, || and [L^] . Since 
we have analyzed the same data with improved 
theoretical input jjj, it is worthwhile to under- 
stand the origin of the numerical difference. 

ALEPH makes a global fit to the five measured 
moments, including the last two which are unre- 
liable (f00% theoretical errors). In view of the 
asymptotic behaviour of A^ ; (a T ), they truncate 
this perturbative series at 0(a T ), neglecting the 

known and positive O(a^) and O(a^) contribu- 
te) 

tions. Thus, they use a smaller value of A y kl '(a T ) 
and, therefore, get a larger result for m s (the first 
value above) because the sensitivity to this pa- 
rameter is through the product m^(M^)A^ (a T ). 
Since they put rather conservative errors, their 
result is nevertheless consistent with ours. 

ALEPH has made a second analysis subtract- 
ing the J = L contribution. Unfortunately, 
only the pion and kaon contributions are known. 
Using the positivity of the longitudinal spectral 
functions, this pole contributions provide lower 
bounds on Imn^(s) and Indies), which trans- 
late into lower limits on the corresponding J = L 
contribution to SRf. Subtracting this contribu- 
tion, one gets upper bounds on SRf L+T jjj which 
imply m s (M?) < 202 MeV §. 

However, besides subtracting the pion and kaon 
poles, ALEPH makes a tiny ad-hoc correction to 
account for the remaining unknown J = L con- 
tribution, and quotes the resulting number as a 
m s (M%) determination [the second value above]. 
Since they add a generous uncertainty, their num- 
ber does not disagree with ours. However, it is 
actually an upper bound on m s (M%) and not a 
determination of this parameter. 
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